Perturbative Expansion around the Gaussian Effective Potential 

of the Fermion Field Theory 

Geon Hyoung Lee* Tack Hwi Lee^ and Jae Hyung Yee* 
Department of Physics and Institute for Mathematical Sciences 

Yonsei University 
Seoul 120-749, Korea 

Abstract 

We have extended the perturbative expansion method around the Gaus- 
sian effective action to the fermionic field theory, by taking the 2-dimensional 
Gross-Neveu model as an example. We have computed both the zero tem- 
perature and the finite temperature effective potentials of the Gross-Neveu 
model up to the first perturbative correction terms, and have found that the 
critical temperature, at which dynamically broken symmetry is restored, is 
significantly improved for the small value of the flavour number. 



*email: ghlee@theory.yonsei.ac.kr 
^ email: thlee@theory.yonsei.ac.kr 
^ email: jhyee@phya.yonsei.ac.kr 



1 



I. INTRODUCTION 



One of the most important problems in quantum field theory is to develop efficient meth- 
ods of non-perturbative approximation. The variational Gaussian approximation method 
has provided an efficient and convenient device in obtaining non-perturbative information 
from various quantum field theories ffl]. The variational method, however, has a serious 
shortcoming that it does not provide a systematic technique to compute the correction 
terms to the approximation. Efforts to establish systematic method to improve the varia- 
tional approximation have been made by the authors of Ref. [H ||. 

Recently a systematic method of perturbative expansion around the Gaussian effective 
action(GEA) [|] has been developed based on the background field method ||. This method 
provides an efficient device to compute the generating functionals for the one-particle- 
irreducible Green's functions in perturbation series, whose zeroth order term is the GEA. For 
the effective potentials of time-independent systems, the result of this method is the same 
as that of the variational perturbation theory developed by Cea an Tedesco ||. It has been 
shown, for the quantum mechanical anharmonic oscillator, that the perturbative correction 
greatly improves the Gaussian approximation even at the first non-trivial correction level 

E- 

It is the purpose of this paper to extend the perturbative expansion method around 
GEA to the case of fermionic field theories. In the next section a brief review on background 
field method is given. We then develop the perturbative expansion method around GEA for 
the two-dimensional Gross-Neveu model. In section we evaluate the finite temperature 
effective potential for the Gross-Neveu model Q, and show how the perturbative correction 
improves the critical temperature, at which dynamically broken symmetry is restored, from 
that of the Gaussian approximation. We conclude with some discussions in the last section. 
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II. PERTURB ATIVE EXPANSION AROUND THE GAUSSIAN EFFECTIVE 
ACTION OF THE GROSS-NEVEU MODEL 



To give a brief review on the background field method ||, we start from the action 

S[<f>]= I ' d d xC[cj>{x),d,cj>{x)l (1) 



in the d-dimensional space-time, where (f) can be a bosonic or fermionic field variable. The 
generating functional for Green's functions is defined by 

< + |0- > J = e iW = J V<pe iS[ip]+iJip , (2) 

where W[J] is the generating functional for the connected Green's functions, J is the external 
source, and the integral convension, J<p = J d xJ(x)(f>(x), is used in the exponent. The 
vacuum expectation value of the field operator in the presence of external source is defined 
by 

and the effective action is defined by the Legendre transformation, 

r[p] = w[J\-<pj. (4) 

The functional derivative of T[<p\ with respect to if gives, 

= -j, (5) 

which is of the same form as the classical equation of motion. 

We now introduce a new action S[(f) + B} obtained by shifting the field by a background 
field B. This new action defines a new effective action T[<f, B], where (p is vacuum expecta- 
tion value of (p field in the presence of the background field B. One can then show that the 
effective action(^) can be represented as, 

r[<p] = r[<p = o,B]\ B= t. (6) 



In other words, the effective action T[ip] can be obtained by summing all the one-y5-particle- 
irreducible diagrams with no external (p lines. This greatly simplifies the perturbative com- 
putation of the effective action. 

It has been shown in ref. that one can rearrange the diagrams in the effective action 
in such a way that the propagator used in the perturbative expansion becomes that of 
the Gaussian approximation, and the zeroth order term of the effective action consists of 
the GEA. 

To extend this method to the case of fermionic field theories, we consider the 2- 
dimensional Gross-Neveu model || described by the Langrangian density, 

C = ipaifiipa - g4> a ip a a - ^cr 2 , a = 1, 2, 3...N, (7) 

where ip a is 2-dimensional Dirac field with a denoting the flavour component. The generating 
functional for Green's functions is then given by, 

Zfo, fj,J) = J Vi>V^Vae l I d2xCo+i ^ +i ^ +l(7j . (8) 

We are interested in the effective action as a functional of the vacuum expetation value of 
the composite field operator, < o >= —g < ipip >, and thus we shift only the composite 
operator; a — > o + B. Then the generating functional for the Green's functions in the 
presence of the background field B is given by 

Z{r],fj, J, B] = J VipVipVae 1 / d2x ^^-h( a+B ) 2 ~9^( a+B ) +i ^ +i ^ T > +iaJ 

where 

K~ l = -(p + igB) 

Kb 1 = -i. (10) 
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The first three factors of the last line of (^|) represent the one-loop effective action, and the 
remaining factor, upon setting a = ^jjlogZ[r],fj, J, B] = 0, gives the higher loop contribu- 
tions to the effective action. Note that the fermion propagator in Eq.([T0|) already contains 
the interaction effect through the background field B, which is the reason why the procedure 
of computing higher loop contributions is simplified. 

In order to rearrange the generating functional so that the functional derivative terms 
in Eq.(|9]) represent the perturbative expansion around GEA, we follow the procedure of ref. 
Q. To do this, we consider the following relation, 

J^W x ^ GT,+yGBJ = [GZ + (^)-(^KG*^^ JOflJ - (li) 

where the repeated index x implies the integration over x, and G and Gb are arbitrary 
two-point Green's functions for fermionic and bosonic fields, respectively. Note that this 
type of functional derivative appears in the expansion of the functional derivative factor in 
Eq.(H). Eq.flTTl) contains a diagram where an internal line coming out of a point goes back 
to the same point (cactus-type diagram) which constitues GEA, i.e., the first term in the 
right-hand side of Eq. (|TT|) . 

To extract such cactus-type diagrams out of the perturbative expansion, we define the 
primed functional derivative as, 

(J_J_J_)'= + G iBa! — . (12) 

5r]% 5fj% 5J X 5r/ x 5fj% 5J X aa 5J X 



Then the primed derivative acting on the generating functional becomes, 
,555 



)' e vGvHJG B J = (fjG^G^KGBJfe^ 1 ^^ ', (13) 



5r] a x 5r] x 5J x ' 

which does not contain any cactus-type diagrams. In order to express the generating func- 
tional (||) in terms of the primed derivative, we need to find Green's functions G and Gb 
which satisfy, 

e -9-^-^-h e nKr,+ \ J K B J 



By using the definition of the primed derivative fll2|) , we easily find G,Gb,N and A that 
satisfy 0: 

AT = 1 

A = -gG aa 

G = K 

G B = K B , (15) 



where the new Green's functions G and Gb are the same as those of Eq.(lO). The reason 
why the Green's functions G and Gb are simple in this case is that we have introduced 
the composite operator a, which turns the four fermion interaction into the three particle 



interaction. From Eqs.(||) and(|T4l) one finally obtains the effective action: 



e * r [ B l = Z[r), f], J, B] \ a= 



a 



det{K- x )det{-iK%)e-^ B ' 1 I[B], (16) 



where 



I[B] = e-^x^e-^^&'e^+yxBJ |- =0 . (17) 

Since the functional I[B] does not contain any cactus- type diagrams, the coeffient of I[B] 
in Eq.(|l6D gives rise to the GEA as in the case of scalar 4 theory [|J. Since I[B] can be 
expanded as a power series in the coupling constant g, we have the perturbative expansion 
of the effective action around GEA. 

The linear term in the exponent of Eq. (|17D generates tadpole diagrams, which does not 
contribute to the effetive action ||. Note that Eq.(fP7|) has the same struture as the higher 
order contribution part of (|]) except that Eq.(|TT|) involves only the primed derivative. We 
can therefore compute the perturbative correction terms to GEA using the same procedure 
as the conventional background field method, by using the Feynman rule (in momemtum 
space) , 
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propagators : K = , m = gB, 

7 ■ p + m 

K B = i, 
vertex : — g, 

/d?p 
— -. (18) 

Thus the perturbative correction, I[B], to GEA consists of one a-particle-irreducible 
bubble diagrams with no external lines and without cactus-type diagrams, as in the case of 
4 theory [|]. 

Up to the first non-trivial contribution from I [B], the effective action for the Gross-Neveu 
model can easily be shown to be, 

T = -%Tr ln^- 1 ) - X -B 2 + tg 2 (K B ) xy (K ca ) xy (K ac ) xy , (19) 

where the first two terms are the GEA and the last term is the first order perturbative 
correction to GEA. 

The effective potential is defined by 

V«M - -fjg, (20) 

where B is the space-time independent background field. Thus the effective potential of the 
Gross-Neveu model, up to the first order perturbative correction term, is given by 

V eff (m) = V G (m) + V P (m), (21) 

where Vq is the Gaussian effective potential, 

and Vp is the perturbative correction, 

V P = ggM/^ * ] 2 = ^m 2 (ln^) 2 , (23) 
2 y 7 2tt ^jFTm 21 87r 2 A 2 K ' 

with A being the ultraviolet momentum cutoff. 



To extract physical information from the effective potentials Vg or V e ff, we need to renor- 
malize them. We can renormalize the effective potentials by requiring the renormalization 
conditions, 

— I " - (24) 
go = g R o R , (25) 

where mo represents a renormalization point. 

For the Gaussian effective potential, Vq, the renormalization condition becomes 

£|U*=^ = V£(ln$ + 2). (26) 
dm z g R g l 2tt A z 

Substituting Eq.(p6|) into (|22|), one obtains the renormalized Gaussian effective potential, 

1 m 2 Nm 2 m 2 . , . 

Va = 2A + ^ (Xn <-^ (27) 

which is equivalent to that of the large- N approximation |J. 

For the perturbatively corrected Gaussian effective potential, V e ff of Eq . (|2l|) , the renor- 
malization condition (j25l) becomes 

*? l—= S = ? + ^ T$ + 2 » + A? + 6 ln + *1- < 28 > 

Requiring that the renormalized coupling constant g R be finite, we find the condition, 



? = <-l±^-^)l^,. (29) 



for which the effective potential can be made finite. The case of the negative sign in Eq. 
leads to unphysical theory. Taking the case of positive sign in Eq.(^), we finally obtain the 
renormalized effective potential, 



1 m 2 , T / 4 m 2 m 2 , . 

2 g R V N 4-7T ' mo 

which clearly shows that the higher order corrections to the Gaussian or the large N ap- 
proximation. This shows that our method is meanigful when N is lager than 4 for the 
Gross- Neveu model. 
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III. EFFECTIVE POTENTIAL AT FINITE TEMPERATURE 



To illustrate how our method improves the Gaussian approximation, we compute the 
first non-trivial perturbative correction term to the Gaussian effective potential at finite 
temperature and evaluate the critical temperature, by using the imaginary time formulation 
of the finite temperature quantum field theory. The generating functional for the thermal 
Green's functions is given by, 

Z p [J\ = N [ V(pe J p dtJ dx c(^) +v j^ ( 31 ) 

J periodic 

where r is the imaginary time with period j3, the inverse temperature, (p = and the 
functional integration is performed only for the periodic field (p. Since the finite temperature 
generating functional is of the same form as the zero-temperature generating functional 
(§) except that it is defined in Euclidean space with periodic boundery condition, we can 
evaluate the perturbative expansion by following the same procedure as in the last section, 
except that the Feynman rule is now modified by, 

propagator : , p„ = (p = iu^pt), u n = (2n + l)7r/? _1 

7 • p — m 

i r dp 

loop integration : - ^ / — . (32) 

' i— — oo 

One can then find the finite temperature effective potential, up to the first perturbative 
correction, 

- V + J Nl + + ^[/ - ^)P. (33, 

where = ^ s /p ^ +~g J o yI . 

To renormalize the effective potential (|33D , it is convenient to seperate the zero- 
temperature and the finite temperature parts of the effective potential: 

Vf ff = V + V p , (34) 

where 



v °=? - 2N J£2 + —y£j - (35) 

V S = -2N- R I + gVN[j ^-}F(p,m) + ^-[F(p,m)]\ (36) 
with m = ga and m) is denned by 

F(J3,m)= I ^1—^—. (37) 



We note that the zero temperature part of the effective potential fl35|) is the same as ( |34"D of 
the last section, and we need only to renormalize this part of the effective potential. Thus 
we have the renormalized effective potential at finite temperature, 



where m = gRfjR. Since F(j3, m) is finite, the renormalization condition ( p9|) reduces the 
finite temperature part of the effective potential, Vp, to 



Vp = -2Ari / g ln(l + e~^) |]m 2 F(/5, m). (39) 

The symmetry that is broken at zero temperature is restored as temperature increases 
beyond the critical temperature. It is well known that the critical temperature should be zero 
in 2-dimensional space-time 0, while in the large-N || and the Gaussian approximations 
of the Gross-Neveu model imply the non-vanishing critical temperature. To see how the 
perturbatively improved Gaussian approximation improves this result, we now evaluate the 
critical temperature from the effective potential (|39l ) . 

Dynamical symmetry breaking in the Gross-Neveu model is manifested by the fact that 
the minimum of the zero temperature effective potential occurs at non-vanishing value of 
the composite field, <jr = cx m , which breaks the symmetry of the classical potential, i.e., the 
symmetry under a — > —a. a m is determined by, 



^■^f^}-^- 1 (40) 
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As the temperature increases the value of a m (P), at which Vfj^ is minimized, decreases. At 
the critical temperature, T c = 4-, cr m {(3) vanishes, which implies the symmetry restoration. 
Thus the critical temperature is determined by 

dV 

— \a R =a m (p c )=0= 0, (41) 

which, together with Eqs. (|38|) and fl39|), gives the value of the critical temperature, 



T c = 0.57g R a m . (42) 
Due to Eq.(f40|), this can be written as, 

7T 1 

T c = 0.57g R a exp[l T ], (43) 

where o"o represents the renormalization point. The critical temperature in the Gaussian 
approximation can similarly be obtained : 



rpgaussian q ^Yg {j9 auss,lan 

7T 1 

7 r n 



Q.57g R a exp[l - —). (44) 



Although the critical temperature T c from the perturbatively corrected Gaussian approx- 
imation does not vanish, it is smaller than that of the Gaussian approximation. This shows 
that the perturbatively improved Gaussian approximation significantly improves the critical 
temperature for small N, the flavour number of the fermion field even at the first order 
correction level. 



IV. DISCUSSION 

We have extended the method of the perturbative expansion around GEA developed 
in ref. [|J to the fermionic field theory, taking the Gross-Neveu model as an example. This 
method is based on the observation that the Gaussian effective action consists of cactus-type 
diagrams, which is extracted out of the functional derivative part of the effective action, i.e., 
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the last two factors of the last line of Eq.(|9|), by introducing the primed functional derivative 
defined in Eq. (|T2"D . This procedure effectively rearranges the diagrams in such a way that 
the zeroth order term of the effective action consists of the GEA. Due to the introduction 
of the composite field, a = —gtptp, as an order parameter, the expansion works only for the 
flavour number N larger than 4. 

In the last section the finite temperature effective potential of the Gross-Neveu model 
is obtained up to the first non-trivial perturbative correction terms, and the critical tem- 
perature, at which the dynamically broken symmetry is restored, is evaluated. The result 
shows a significant improvement of the critical temperature compared to the Gaussian and 
the lagre-iV results for small value of the flavour number. 
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